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I. INTRODUCTION 



The meson light-cone distribution amplitude (DA) is of great importance to understand 
the key properties of mesons in exclusive hardronic processes, since it provides essential 
information on the non-perturbative structure of a meson in the processes 0, 0, S 0]- An 
experimental aspect of the pion DA can be seen in the recent measurements of 7r-photon 
transition form factor by the CLEO experiment (if . The analysis of these data in Ref . 
has shown that neither two-humped DA for the pion predicted by Chernyak and Zhitnitsky 
(CZ) 0| nor the asymptotic one are favored at the 2a level of accuracy. These findings 
were independently confirmed by Bakulev et al. in Ref. jHEJ. On the theoretical side, the 
pion DA has been investigated in various approaches: For example, in the QCD sum rules 
(QCDSR) H E3 Q 0, El , in lattice QCD fLQCD) [3 Gil, in the chiralquark model 



from the instanton vacuum 0, HEl 0, H 0,0, in the NJL models 0, 0, £gj, m 
chiral perturbation theory (yPT) \2l\ and so on. A comparison of the results from those 
models with the CLEO data was lucidly made in Ref. 0], all the results being evaluated at 
/i 2 = 5.76 GeV 2 (known as the Schmedding-Yakovlev (SY) scale after NLO evolution. 

While the pion DA has been extensively studied for well over decades, the kaon one has 
attracted attention rather recently, since it is of deep relevance to the exclusive 5-meson 
decays to light mesons. The 5-meson decays are now being under full investigation at the 
BaBar and Belle experiments. These experiments will soon shed light on the pattern of the 
CP-violation as well as of flavor SU(3)-symmetry breaking. In this context, there has been 
some amount of theoretical works on the DA in the QCDSR 3(1 U, 3^, 33l 34| . Though CZ 
derived kaon DA, originally in the QCDSR 0, 0], their results are known to suffer from 
the sign mistake. 

It is known that spontaneous chiral symmetry breaking is well realized in the instanton 
vacuum via quark zero modes. Thus, it may provide a good framework to study the structure 
of light mesons. Moreover, there are only two parameters in this approach: The average 
instanton size p ~ | fm and average inter-instanton distance R ~ 1 fm. In particular, the 
normalization scale of the present approach is naturally identified by the average size of 
instantons and is approximately equal to p _1 ~ 600 MeV. The values of the p and R were 
estimated phenomenologically by Ref. 0] as well as theoretically by Refs. 0, 0] . These 
values were recently confirmed in various lattice simulations of the QCD vacuum 0, 41 



12]. Very recent lattice calculations of the quark propagator 0, 0| are in a remarkable 



agreement with that of Ref. [38]. The instanton vacuum model was later extended by 
introducing the current-quark masses 0, Ull 0|. It was assumed in the model that the 
large N c expansion is the reasonable one and the results were obtained in the leading order in 
this expansion. In the present work, we want to investigate the leading-twist pion and kaon 
DAs within the framework of the nonlocal chiral quark model (xQMj from the instanton 
vacuum, with SU(3)-symmetry breaking effects taken into account |45|, |46|, |4jj . We basically 
follow the formalism in Ref. |24| apart from flavor SU(3)-symmetry breaking effects. 

Although the derivation of the effective action with flavor SU(3)-symmetry breaking 
was rather complicated, the final result is summarized in the denominator of the quark 
propagator with explicit flavor SU(3)-symmetry breaking in a very simple form, ft + i[mf + 
Mf(k)}. The only difference is the explicit appearance of the current-quark mass in the 
denominator when it is compared to the propagator in the chiral limit (m s = m u = 0). 
We note here that Mf(k) is the dynamical quark mass depending on the quark momentum 
and current-quark mass. However, instead of using the non-local quark mass derived from 
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the instanton vacuum, we employ in the present work a simple-pole parameterization of the 
momentum-dependent quark mass Mf(k), since it shows a similar behavior to the original 
one and it can be easily continued analytically to Minkowski space. 

The xQM from the instanton vacuum has one great virtue: The normalization point is 
naturally given by the average instanton size whose value is determined by the saddle-point 
equation. However, we will fix the scale by the normalization condition for the DAs in 
the present work in such a way that we can reproduce the pion and kaon decay constants 
simultaneously. Employing the dipole-type form factor, we get the cut-off mass A to be 
around 1.2 GeV which plays a role of the scale of the present method. However, since the 
evolution does not depend sensitively on the scale, we will take our scale to be 1 GeV which 
is not far from the original scale [17[. From this value, we can evolve the DAs obtained in 
this approach to the scale of the CLEO experiment, A ~ 2.4 GeV (the SY scale [6j) so that 
we can compare the results with the empirical data. 

We will provide in the present work the following numerical results: The DAs, the Gegen- 
bauer moments, the moments of the amplitudes, the transverse momentum distributions and 
their ratios for the pion and kaon. We observe the symmetric and asymmetric shapes for the 
pion and kaon, respectively, as expected from the effects of explicit flavor SU(3)-symmetry 
breaking. However, when we take into account the current-quark mass corrections to the 
dynamical quark mass, the kaon DA becomes less asymmetric. We also consider the one- 
loop QCD evolution of the DA so that we may compare our results with those from the 
other models. Our results are in good agreement with others qualitatively. The moments of 
the distribution amplitudes, (£ m ), are also given for the pion and kaon. Finally, we consider 
the transverse momentum distributions of the full light-cone wave function, Using 
this, we test the factorization hypothesis for the kaon. 

This paper is organized as follows: In Section II, we briefly explain the general formalism 
to derive the DAs within the framework of the nonlocal xQM from the instanton vacuum. 
The numerical results are given in Section III with discussions. The last Section is devoted 
to summary and conclusion. 



II. GENERAL FORMALISM 



In the present Section, we briefly explain how to derive the pion and kaon light-cone DAs 
within the framework of the nonlocal xQM from the instanton vacuum. The leading-twist 
(twist-two) light-cone DA for the pion and kaon is defined as follows 0]: 



00 dz 
— < 
-oo 2ir 



- i{2u - 1)z - p (0\q f (z)fi lb exp 



ig / dz"*AJz') 



q g (-z)\<f>(P)), 



where \<fr(P)) is the state vector for a pesudoscalar meson moving with on mass-shell momen- 
tum P in the light-cone frame, ip and stand for the quark and gluon fields, respectively. 
u and z indicate the longitudinal momentum fraction and spatial separation between the 
quarks, is the light-like vector satisfying the condition of n 2 = 0. The exponential term 
is called the Willson line or gauge connection which guarantees the gauge invariance of the 
nonlocal quark bilinear operator. However, by virtue of the light-cone gauge, we set this 
Wilson line to be unity. Therefore, we need not consider any gluon contribution for leading- 
twist pion and kaon DAs. The pseudoscalar (PS) meson decay constant is introduced as 
the normalization constant, which will be determined with the cut-off mass A in the model. 
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The DA must fulfill the following normalization condition: 



1. 



(2) 



We will fix the cut-off mass A which is approximately identified as the scale of the present 
approach, using Eq. (J2J). In order to calculate the pion and kaon DAs, we first start with 
the partition function of the nonlocal xQM in Euclidean space: 



Z I Vi[)V^V4> a exp j d A x[^f{x){ifi + im f )ipj 
d A kdH 



x 



+ i 



(2tt)> 



e ^-D-^ Mf (k)M g (l)^(k) (U^) fg ^(l) 



(3) 



This partition function can be constructed by considering arbitrary Nf, large N c and the 
saddle-point equation based on the instanton induced 2iVj-quark interaction. The subscripts 
/ and g denote the quark-flavor indices. In the present work, we have (/, g, <fi) = (s, u, K + ) 
and (d,u,ir + ). a represents the color index. Note that the quark and anti-quark interact 
each other nonlocally via the nonlinear PS-background field of <p a . The dynamical quark 
mass Mf(k) plays a role of the quark and PS-meson coupling strength, which arises originally 
from the Fourier transform of the quark zero-mode solutions. Moreover, this momentum- 
dependent quark mass can be regarded as a natural and intrinsic UV regulator for loops in 
the model. The nonlinear background PS-meson field is defined as follows: 



exp(* 75 a (x)A7^) 



(4) 



where A a is the flavor SU(3) Gell-Mann matrices, the quark propagator can be easily ob- 
tained from the effective chiral action given in Eq. (jSJ). One of the authors (M.M.) suggested 
the effective chiral action in terms of the multi-quark interaction with explicit flavor SU(3)- 
symmetry breaking ji^]. This action is called the modified improved action (MIA), which 
can be distinguished from the usual effective action derived by Diakonov and Petrov I3T 
The quark propagator from the MIA is obtained as follows: 



1 



+ i[m f + M f (k)} 



(5) 



Note that the denominator contains the current-quark mass rrif explicitly. Moreover, the 
dynamical quark mass Mf(k) also depends on the m/. We parameterize in the present work 
the dynamical quark mass as follows: 



M f (k) = M F 2 (k 2 )f(mj 



M n 



nA 2 



(nA 2 + k 2 ) 



2)i 



1 + 



d 2 



mj_ 
d 



(6) 



This parameterization of a simple-pole type resembles the behavior of the Fourier transform 
of the quark zero modes. M is the constituent quark mass being set to be 350 MeV. 



24l |. The form of the 



We consider the power n in Eq. (|6J) free parameter 

correction factor /(m/) in Eq. © was suggested by Ref. [45| via the saddle-point equation 
of the effective action. The explicit expression of f{rnf) was derived by Pobylitsa: If one 
expands quark propagator in the presence of the instanton effect and resums quark loops 
which are suppressed in the large N c expansion, one gets the correction factor fijrif) as 
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given in Eq. (JBJ) [48j]. For convenience, we denote the effective action without the correction 
factor by MIAl and that with it by MIA 2. 

It is worth mentioning that there is a caveat in the present approach. The effective 
low-energy QCD partition function of Eq. (jSJ) is well defined in Euclidean space. However, 
we have to continue analytically to Minkowski space, in order to study the light-cone DAs. 
Though we have no theoretically firm grounds for such analytic continuation, we want to 
adopt a more practical stance on it. Since the nonlocal xQM from the instanton vacuum was 
used for investigating the light-cone DAs with this analytic continuation successfully [l7| . 
we will proceed to study the DAs in Minkowski space very carefully. 

The DAs in the present approach is obtained as follows: 



d A k 



5(2uP ■h-2k-n) 



x tr 



M f (k) 



M f (k-P) 
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r d 4 k 
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(?) 



where the functions F, G and D are given as follows: 

F f = (k 2 - nA 2 ) 3n , 
Gf 
Df 



(k - nA 2 ) n [m f (k 2 - A 2 ) 2n + rj\ 



nA 



(k - nA 2 Y n (k 
[(k - P) 2 
[(k - P) 2 
[(k - P) 2 



m 



- 2m f rj(k 2 - nA 2 ) 2n - r] 2 + ie 



2l3n 



nA 2 ] n [-m g [(k - P) 2 - nA 2 } 211 - V ], 



nA 2 ] 4n [(A; - Pf 



2m g r][(k -P) 2 - nA 2 } 211 - V 2 + ie 



with r] = M A 4n . We evaluate Eq. © to order 0( 



m 



using the relations P 2 
P£ ~ P + P~ = m 2 K in the light-cone frame. We set = 140 MeV and mx 
for numerical input. 



: p+p- - 
495 MeV 



III. NUMERICAL RESULTS 

First, we show briefly how to fix the cut-off mass A. We set A by using the normalization 
condition given in Eq. (J2J) in such a way that it reproduces the empirical tt and K meson 
decay constants simultaneously. We first consider the case of n = 1. In this case, we obtain 
A = 1.2 GeV with which the pion decay constant F n = 96.77 MeV (93 MeV) and the kaon 
decay constant Fx = 110.21 MeV (113 MeV) are well reproduced qualitatively in comparison 
with the experimental data put in the parentheses. We use the value of A = 1.2 GeV for all 
cases of n and list the results of the calculated pion and kaon decay constants in Table HI As 
shown in Table HI we find that the values of the pion and kaon decay constants decrease as n 
increases. In the case of the MIAl, they are qualitatively in a good agreement with the data, 
while for the MIA2 the kaon decay constants are underestimated by about 20 %. It can be 
understood as follows: The current-quark mass corrections reduce Mf(k,rrif = 0) by about 
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25% for the strange quark whereas they lessen it by about 2% for the up and down quarks. 
We will see later that the kaon distribution amplitudes is also reduced similarly with the 
MIA2 used. Though we are able to fit the kaon decay constant for the MIA2 by changing 





F v 


F K 


n 


12 3 


1 


2 3 


MIA1 


96.77 93.66 92.83 


110.21 


107.22 106.36 


MIA2 


95.13 92.45 91.76 


98.51 


96.47 96.07 



TABLE I: The results of the pion decay constant F n and the kaon decay constant Fk in units of 
MeV with A = 1.2 GeV. 

the cut-off mass, we will not try to do that, since we want to produce all results with one 
cut-off mass. Moreover, note that the meson-loop corrections (~ 1/AQ are not taken into 
account in the present work, which is known to be of great significance to reproduce the 
proper value of the kaon decay constant as shown in as well as in the nonlocal xQM 
calculations [49l l50t l5l|. However, we do not take into account any 1/N C correction in the 
present work, since it is beyond the aim of the present work. Thus, we still have room to 
improve the result of the kaon decay constant in the future |Hl| . including the meson-loop 
corrections. 

The presence of the nonlocal interaction between quarks and Nambu-Goldstone bosons, 
which arises from the momentum-dependent quark mass, breaks the gauge invariance, so 
that the currents are not conserved. The conserved currents in Euclidean space with the 
nonlocal interactions can be derived by gauging the partition function [H^. The correct 
expression of the pion decay constant / 2 can be derived by using the modified axial-vector 
current in the following matrix elements: 

7T 6 (P)} = if^P^S ", (8) 

which indicates that the Ward-Takahashi identity of PCAC is well satisfied with the modified 

conserved axial- vector current. If we use the usual form of the axial- vector current Af, = 

_ p 

ip r Yti'j5^ a "ip, we would end up with the Pagels-Stokar (PS) expression for / 2 PS : 
r d^k M 2 — -MM'k 

/£ps = 4N c J ( 27r )4 (P + M 2 ) 2 (in Euclidean space) - (9) 

which gives smaller results than the correct expression by approximately 20 %. Thus, one 
has to consider the modified conserved currents in Eq. (j7J). However, if we use the /^pg for 
the normalization of the pion and kaon DAs for convenience, we need not take into account 
the current conservation in Eq. (j7J), since the results are almost the same jH^. 

In Fig. ^ we depict the pion and kaon DAs, varying the power n with two different 
modified improved actions, MIAl and MIA2. The asymptotic form of the distribution 
amplitude $A sy m = 6u(l — u) is plotted in each panel of the figure for comparison. The 
results of the pion DA are all symmetric under the inversion u —>■ 1 — u due to the fact that 
the mass of the quarks inside the pion are negligibly sm all ( m,,/A ~ m^/A ~ 0). Thus, 
the present results are almost the same as those in Ref . [2J| (only 2 % difference) . Note 
that the current-quark masses for all flavors are taken into account in the present work. All 
curves are rather flat in the vicinity of the centered area (0.25 ^$ u ^ 0.75), compared to 
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MIA1 



MIA2 




1.00 



1.00 



FIG. 1: DAs for the pion (upper two panels) and kaon (lower two panels) with different powers n. 
In the right column, we draw the distribution amplitudes for the MIAl, while we depict in the left 
column those for the MIA2. We also draw the asymptotic distribution amplitude in each panel for 
comparison. 



the asymptotic one as already observed in Ref . [13[ . The pion DAs vanish at the end-points 
u = and 1, which indicates that it is impossible for a quark in the pion to have zero or 
full momentum of P. In the present approach, this behavior is caused mainly by the term 



with yj Mf(k)M g (k) in the DA in Eq. (J7|). It infers that the nonlocal interaction between 
quarks and pseudoscalar mesons is essential in describing the correct end-point behavior of 
the DAs. 

Interestingly, as the power n increases, the results of the pion DA behave differently at 
the end-points, u = and 1. The results are getting suppressed as the n gets larger, so 
that the end-point behavior of the DA changes from a convex form to a concave one. This 
distinctive behavior is of great importance for comparison with the empirical data 0, 1 . 



See Refs. for some relevant discussion and for the role of the dynamical quark mass for the 

end-point behavior. 
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The kaon DAs are drawn in the lower panels of Fig. Q]for the MIAl and MIA2. The overall 
shapes and behavior at the end-points are very similar to those of the pion ones. However, 
as expected from explicit flavor SU(3)-symmetry breaking, the results become asymmetric, 
since the quarks inside the kaon carry different fractions of the total momentum P due 
to the mass difference between the strange quark (m s = 150 MeV) and the up and down 
quarks {m U( i = 5 MeV). We find that there is a noticeable difference between the cases of 
the MIAl and MIA2. As for the MIAl, the results look more asymmetric, i.e. the right 
side is larger than the left side whereas those for the MIA2 are almost symmetric. This 
can be understood by recalling the fact that the current-quark mass correction factor f(rrif) 
reduces the dynamical quark mass by ~ 25%, so that the mass difference between the light 
quark and strange quark gets smaller: 

m u + M u (k = 0) =ori ^ m u + M u (k = 0) = Q ^ 
m s + M s (k = 0) MIA1 ' m s + M s (k = 0) MIA2 

Thus, the DAs for the MIA2 turn out to be nearly symmetric. We note that the results of 
the kaon DA for the MIAl are rather similar to those calculated without the tensor term 
~ kfjc* in Ref. @. 

We are now in a position to study the Gegenbauer moments for the pion and kaon DAs. 
The Gegenbauer moments are very useful in analyzing the DAs. It can be extracted from 
the experiments or derived from theoretical models. In principle, the Gegenbauer moments 
indicate how much the DAs deviate from the asymptotic one. They are defined by the 
coefficient of the expansion of the PS-meson DA in terms of the Gegenbauer polynomials: 

oo 

$(u) = $ Asym W E a-Cm /2 (£) 

m=l 

= 6u(l - u)[a C 3 /2 (0 + a x Cl' 2 {0 + a 2 C 2 3/2 (0 + ■•■], (11) 

where aoC^ 2 (^) = 1 and £ = 2u — l. The coefficients a m are called the Gegenbauer moments 
and can be calculated by using the orthogonal condition for the Gegenbauer polynomials: 

1 (i - ^cT (0c£?(0 = r j^j- s mm , (12) 

-i ml{m + S/2) 

In Table |Hj we list the calculated Gegenbauer moments for the pion DAs with n varied for 
the MIAl and MIA2. Due to isospin symmetry, the odd Gegenbauer moments of the pion 
DA become zero. Therefore, we only need to consider the even moments. In spite of the fact 
that u- and <i-quarks have slightly different masses (m u ~ 5 MeV and ~ 10 MeV), the 
difference can be ignored in comparison to the scale of the model, A ~ 1 GeV. We observe 
the second Gegenbauer moments are gradually decrease as the n increases without sign 
changes for both of the MIAl and MIA2. On the contrary, the fourth and sixth Gegenbauer 
moments turn negative as the n grows. It is deeply related to the fact that the end-point 
behavior of the pion DAs is changed as the n gets larger. These sign changes for the fourth 
and sixth Gegenbauer moments govern the shapes of the pion DAs around the end-points. 
We also list the results from other theoretical calculations, mainly from the QCDSR and 
lattice QCD with specific scales denoted in the brackets. Since our scale is about 1 GeV, in 



order to compare the present results with those from Refs. p, |26|, |5JJ at the SY scale, we 



need to evolve the present results to the corresponding scales. The anomalous dimension 
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n 


(2<2 ^4 


1, MIA1 

C\ TV /TT A 1 

2, JV11A1 

3, MIA1 


0.11012 0.01500 0.00153 
0.05330 —0.02831 —0.01597 
0.02903 -0.04603 -0.01897 


1, MIA2 

2, MIA2 

3, MIA2 


0.11399 0.01583 0.00188 
0.05963 -0.02822 -0.01523 
0.03442 -0.04479 -0.01917 


1 

[6j (2.4 GeV) 
I2bl (2.4 GeVJ 
[32] (1.0 GeV) 
[57] (2.24 GeV) 
[fig] (1.0 GeV) 
[59] (1.0 GeV) 


0.56 
0.12 ±0.03 

0.044 ± 0.016 0.023 ± 0.010 

u - ZD -0.09 

0.236(82) 
0.25 ±0.15 

0.44 0.25 



TABLE II: The Gegenbauer moments for the pion distribution amplitudes. The numbers in the 
parentheses stand for the scales of the corresponding works. 



of the Gegenbauer moments for the one-loop renormalization equation can be written as 
follows: 



7 



(0) 



3± 



We have 7 J 0) = 7.11 

moments in two different renormalization scales can be compared by the following relation: 



(m + l)(m + 2) 



m+l i 

- 4 £± 

m'=l 



nv 



7?> 



11.11. 



4 0) 



13.95. 



(o) 

74 



(13) 



16.18, • • •. Thus, the Gegenbauer 



a m (Ai) = a m (A 2 ) 



a(A 2 



7^V(2/3 2 ) 



a m (A 2 ) 



ln[A 2 /A ] 
ln[A x /A ] 



7^/(2/32) 



(14) 



where (3q = 9 and A — 200 MeV. Using Eq. (jl4j) . we can easily compare our results with 

those at 2.4 GeV by multiplying them by 0.72 7 ™ ^ 2/3 °). We see that the value a 2 = 0.09 ~ 0.15 
from Ref. is quite consistent with the present results for the MIAl ranging 0.07 ~ 0.13 
incorporated with the QCD evolution correction, 0.878. 

Now, we consider the kaon Gegenbauer moments. As for the even Gegenbauer moments, 
the general tendency is very similar to those of the pion: The negative sign appears as the 
n grows, which makes the DAs suppressed around the end-points. The results of the second 
Gegenbauer moment a 2 are comparable with those from Refs. |30| . 

The most interesting point can be found in the odd Gegenbauer moments which are finite 
due to the effects of explicit flavor SU(3)-symmetry breaking, as shown in Table ILTT1 As for 
af', we get all positive values for the MIAl whereas they are almost negligible for the MIA2. 
As a result, while the kaon DAs turn out to be asymmetric for the MIAl, they are almost 
symmetric for the MIA2 as shown already in Fig. ^ Ref- [3 i f p rovides the negative value 
with large uncertainty. The positive was suggested in Ref. |32| ranging 0.04 ~ 0.09 which 
is about five times larger than those for the MIA2. When we compare our results with those 
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ra,MIA 


af af af af af af 


1, MIA1 

2, M1A1 

3, MIA1 


0.06865 0.03264 -0.00559 0.00469 0.00315 -0.00031 
0.08828 —0.02634 —0.00657 —0.02820 0.00022 —0.00957 
0.09558 -0.05145 -0.00777 -0.03999 -0.00214 -0.00965 


1, MIA2 

2, MIA2 

3,1V11A2 


-0.00667 0.03348 0.00160 0.00895 -0.00158 0.00088 
-0.00047 -0.01229 0.00010 -0.01928 -0.00248 -0.01017 
U.U0107 — U.031U4 — 0.U0U83 —0.03103 —0.00278 —0.01120 


[31] (1.0 GeV) 
[22] (1.0 GeV) 
[58] (1.0 GeV) 
[60] (1.0 GeV) 


0.10 ±0.12 

n n7 +0.02 n o7 +0.37 
U - U '-0.03 u - z '-0.12 

0.30 ±0.15 

0.06 ± 0.03 



TABLE III: The Gegenbauer moments for the kaon distribution amplitude. The numbers in the 
parentheses stand for the scales of the corresponding works. 



of Ref. jU H E| in which various techniques in the QCDSR were taken into account, 
the MIAl gives qualitative agreement with them for af whereas the MIA2 provides much 
smaller af . The calculated values of af are all much smaller than those estimated from the 
QCDSR. 

In addition to the Gegenbauer moments, we can also define the expectation value of the 
momentum which is called £ -moment, as follows: 

r%= tduc%{u) (15) 







with £ = 2u— 1. The numerical results are listed in Table ITV1 We compute the ^-moments 
upto m = 6 for pion and kaon DAs. As expected from isospin symmetry, the odd £- 
moments for the pion DAs vanish. As for the kaon DAs, we observe all negative odd 
^-moments for the MIAl whereas positive for the MIA2, which is a similar situation in the 
case of the Geg enbauer moments. We also list the ^-moments from other model calculations 
flil 26l . 57, 61, 6^, 6^ for comparison. 



Now we discuss the inverse moment / |28( and may relate to a physical observable in the 
CLEO experiment of 7*7* — > vr° in terms of the pion form factor |5[ : 

du , , . 16 
4« 1-m v ; 



As shown in Ref. [2j| , the value of I turns out to be exactly 3 with the asymptotic distribution 
amplitude. In Table 13 we list the results of / for the MIAl and MIA2 with the n varied. 
Although the results are similar to that of the asymptotic one (1 = 3), they in fact deviate 
from it by approximately 10%. Note that I n and 1% are very similar each other (I K < /„.). 

We now discuss the expectation value of the transverse momentum (k™)- As indicated in 
Refs. jzllill' it contains important information on the meson in hard processes. Moreover, 
it is deeply related to the QCD condensates The moments of $(/ct) can be written as 
follows: 

= T W ^(0\q f HKDTqMP)), (17) 
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n, MIA 




(e) K <ev <ev <ev <ev 


1, MIA1 

2, MIA1 

3, MIA1 


0.23776 0.11244 0.06663 
0.21828 0.09496 0.05206 
0.20995 0.08757 0.04606 


0.04119 0.21119 0.01659 0.09366 0.00902 0.05326 
0.05297 0.19097 0.02145 0.07676 0.01149 0.03985 
0.05735 0.18236 0.02310 0.06980 0.01219 0.03452 


1, MIA2 

2, MIA2 

3, MIA2 


0.23908 0.11341 0.06734 
0.22044 0.09641 0.05307 
0.21180 0.08893 0.04704 


-0.00400 0.21148 -0.00141 0.09430 -0.00076 0.05394 
-0.00028 0.19579 -0.00010 0.08090 -0.00019 0.04309 
0.00064 0.18936 0.00012 0.07540 -0.00015 0.03873 


[15] (2.67 GeV) 
[26] (2.4 GeV) 
[57] (2.24 GeV) 
[61] (1.0 GeV) 
[£2] (1.0 GeV) 
[63] (1.0 GeV) 


0.280i°;°?° 
0.040 ± 0.005 
0.281(28) 
1.37 ±0.20 
0.25 ±0.10 
0.10 ±0.12 





TABLE IV: ^-moments for the pion and kaon DAs. The numbers in the parentheses stand for the 
scales of the corresponding works. 





In 


Ik 


n 


1 2 3 


1 2 3 


MIA1 


3.40471 3.00196 2.88593 


3.12749 2.79872 2.70349 


MIA2 


3.42072 3.02038 2.90240 


3.14755 2.85608 2.77235 




TABLE V: The results 


of the integral I. 



where = — is the covariant derivative in which Z?^ = d +igA tl -\/2, and stands 
for the unit vector in the transverse direction. We note that the kx dependence of the meson 
light-cone wave function \T/ { u, kr ) is derived from the quark propagator (k 2 = k + k_ — h^) as 
shown in Eq. In Refs. |7|,|6J] the (k™) is approximated in terms of quark and quark-gluon 
mixed condensates with the soft pion theorem in the chiral limit: 

2 5 (igqa^G^q) 

(M - 36 — m — ' m 

Hence, there are two different ways to compute {k™ ) in the chiral limit: Firstly, it can be 
computed by using Eq. (fTHj) as done in Ref. 65| within the same scheme as the present work. 



Secondly, it can be directly calculated as follows: 

/*oo rl /*oo 

(k%)= d 2 k T k^${k T )= du d 2 k T k^(k T ,u). (19) 
Jo Jo Jo 

Since the light-cone wave function ^(fcy, m) can be obtained by integrating Eq. (0) over 
k + and k~ , it is possible to calculate (k™) given in Eq. (|19|. In the present work, we only 
consider the cases m = 2 and 4. In Fig. |2l we depict the numerical results for (k 2 ^) 7T) K- As 
for the pion, the transverse momentum is rather equally distributed between the quark and 
the anti-quark inside it, so that (k 2 -,)^ is almost flat within the range 0.25 ^ u < 0.75. On 
the contrary, (fc^)x shows asymmetry, since the light quark and the strange quark inside 
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MIA1 MIA2 




FIG. 2: (A;f.) for the pion (upper two panels) and kaon (lower two panels) for different n. 



the kaon carry different fractions of the transverse momentum due to their mass difference. 
Using the results of (k^(u)), we test now the factorization hypothesis of the light-cone wave 
functions for the pion and kaon: The light-cone wave function can be factorized into the 
longitudinal and transverse contributions in a separable form ^f(u, k?) = (A>r). This 

hypothesis is useful when the meson light-cone wave function is to be integrated over the 
quark four- momentum. However, as discussed in Ref. the hypothesis does not work 
well. It can be understood in the following equation: 

{k 2 T {u)) _ Jd 2 k T kl m ^(k T ,u) fd 2 k T k 2m ${u)<5>{k T ) _ fd 2 k T k 2m §(k T ) 
Wu) ~ Jd 2 k T V(k T ,u) ' Jd 2 k T $(u)$(fc T ) ~ fd 2 k T $(k T ) ' 

(20) 

Thus, if the factorization hypothesis had worked well, then the quantity (k^(u)) / $(w) would 
have been independent of u. In the upper two panels of Fig- El we show the numerical results 
of (k 2 ^ («))/$ (u) for the pion. As shown in Fig. 01 the results are rather flat for the region 
0.25 < u < 0.75 but not constant, which is similar to those of Ref. [25|. In the kaon case, 
the factorization hypothesis is obviously a wrong one. As shown in the lower panel of Fig. 01 
the results are not at all constant. 

Finally, we calculate the averaged expectation values of the transverse momenta with 
respect to u so that we can compare them with those obtained above, using the following 
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MIA1 



MIA2 




FIG. 3: (kfy/Q for the pion (upper two panels) and kaon (lower two panels) for different n. 



equation: 



(k 



,m\ 

T I Averaged 



du dk\ k T ^/(u, kx) 



(21) 



In Table ED we list the results of (k^), (k^) and R = (k^)/(k^) 2 . Interestingly, we get 
rather stable results of the expectation values for n = 2,3. However, we obtain a noticeably 
large result for n = 1 due to the larger values of (k^) for n = 1. As for the kaon, larger 
values of (k^) are yielded, in particular, R ~ 20 with n = 1 for the MIA2. We note that the 
ratio R was estimated to be 5 ~ 9 for the pion DA in Refs. mm, which is similar to the 
present results, especially, in the case of n — 1 (~ 5). However, in general, our results turn 
out to be smaller and less broad than those from other models. We also list those for (k^) n 
from Ref . (26, .64] . The result of Ref . 64] is smaller than ours. 



IV. SUMMARY AND CONCLUSION 

In the present work, we aimed at investigating the leading-twist light-cone distribution 
amplitudes of the pion and kaon, the effects of flavor SU(3)-symmetry breaking being ex- 
plicitly considered from the QCD instanton vacuum. 
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ra,MIA 


(k^n {k^) n {k^ir/ (^t)w 


(k 2 T ) K {k A T ) K {k 4 T ) K /(k 2 T ) 2 K 


1, MIA1 

2, M1A1 

3, MIA1 


0.20671 0.22063 5.16348 
0.17382 0.08618 2.85241 
0.16554 0.07322 2.67164 


0.26937 1.12985 15.57098 
0.21271 0.15226 3.36524 
0.19980 0.11746 2.94253 


1, M1A2 

2, MIA2 

3, MIA2 


n onion n 01 000 c oacoa 

0.20230 0.21882 5.34684 
0.17156 0.08461 2.87481 
0.16302 0.07166 2.69623 


0.22160 0.96318 19.61359 
0.17503 0.12040 3.93006 
0.16391 0.09156 3.40782 


[26J (1 GeV) 
[6j (1 GeV) 


0.18490 
0.10000 





TABLE VI: [GeV 2 ], (4) [GeV 4 ] and (fc|}/(A^) 2 for the pion and kaon DAs for different n 
and MIA. 

We started from the modified improved effective chiral action from the instanton vacuum, 
with two different ways considered: The dynamical quark mass without the current-quark 
mass correction factor (MIAl) and that with it (MIA2). As for the momentum-dependent 
dynamical quark mass, which plays an essential role in describing the meson distribution am- 
plitudes, we parameterized it by a simple-pole type form factor with three different powers, 
n = 1, 2, 3. The cut-off mass was fixed to reproduce pion and kaon decay constants simulta- 
neously in such a way that the normalization condition of the distribution amplitudes may 
be satisfied. 

We first examined the pion and kaon distribution amplitudes in detail. As expected from 
the light current-quark masses {m u d = 5 MeV), we did not find any noticeable difference 
between the pion distribution amplitudes for the MIAl and those for the MIA2. It was 
found that the pion distribution amplitudes turn out to be symmetric for the momentum 
fraction u, while the kaon ones show asymmetric behavior due to flavor SU(3)-symmetry 
breaking effects. However, when the current-quark mass dependence of the dynamical quark 
mass is considered, the kaon distribution amplitudes become less asymmetric. Furthermore, 
we found that as the power n increases, the shape of the distribution amplitudes is getting 
flat. It was also shown that the end-point behavior of the distribution amplitudes varies as 
the power of the form factors is changed. 

In order to see how far the pion and kaon distribution amplitudes deviate from the asymp- 
totic one and how large the flavor SU(3)-symmetry breaking effects are, we calculated the 
Gegenbauer moments of the distribution amplitudes to the sixth order. The odd Gegenbauer 
moments turn out to vanish for the pion distribution amplitudes due to isospin symmetry, 
whereas they remain finite for the kaon one due to the effects of explicit flavor SU(3)- 
symmetry breaking. The fourth and sixth Gegenbauer moments turn negative as the power 
n increases. It implies that the pion distribution amplitude at the end-points {u = 0, 1) is 
changed to be concave, while it becomes flat in the vicinity of u = 0.5, as mentioned above. 
As for the kaon Gegenbauer moments, we observed very similar tendency to those of the 
pion, except for the odd Gegenbauer moments. The current-quark mass corrections to the 
dynamical quark mass made the first Gegenbauer moment nearly negligible. 

We also took into account the QCD evolution via the QCD renormalization group equa- 
tion in order to compare the present results for the Gegenbauer moments with other empirical 
or theoretical estimates at the Schmedding-Yakovlev scale (2.4 GeV). 

Finally, the expectation values of the transverse momentum (k 2 ^ 11 ) were computed by 
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using the distribution amplitudes and light-cone wave functions ty(kT,u). As expected from 
explicit flavor SU(3)-symmetry breaking, the quarks inside the kaon carry the momentum 
fractions in an asymmetric manner, whereas those inside the pion do symmetrically. Using 
these results, we tried to analyze the factorization hypothesis of the total light-cone wave 
function, ^(u, kx) = $(w)$(A>r)- The numerical results of this ratio for the pion look rather 
flat but not constant. Those for the kaon are far from the flat shape. Thus, altogether, 
the factorization hypothesis is invalid for the meson light-cone wave functions. We also 
calculated the momentum fluctuation shown in the ratio of (h?) / '(^f) 2 an d compared it 
with other model calculations. 
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Appendix 

The definition and usage of light-cone coordinate vectors. 

n, = (1, 0, 0, 1), n'^ = (1,0, 0,-1), n-n = 2, 

k + k~ 
k„ = —h^—h'^ + kT^ 

= (0, ki, k 2 , 0), k? = (0, — fci, —k 2 , 0), 

K = ^^(iMi) + ^^(i:OA-i) + (o : k u k 2 ,oi 

*" = ^y^(l,0, 0,-1) + ^±^(1,0, 0,1) + (0,-^,-^,0), 
k + = k + ks = k ■ n, k~ = k — k 3 = k ■ n ', k ■ k = k + k~ — k?, 

d 4 k = -dk + dk-dk 2 T . (22) 



[1] A. V. Efremov and A. V. Radyushkin, Phys. Lett. B 94, 245 (1980). 

[2] G. P. Lepage and S. J. Brodsky, Phys. Lett. B 87, 359 (1979). 

[3] G. P. Lepage and S. J. Brodsky, Phys. Rev. D 22, 2157 (1980). 

[4] V. L. Chernyak and A. R. Zhitnitsky, Phys. Rept. 112, 173 (1984). 

[5] J. Gronberg et al. [CLEO Collaboration], Phys. Rev. D 57, 33 (1998) 

[6] A. Schmedding and O. I. Yakovlev, Phys. Rev. D 62, 116002 (2000). 



15 



[7] V. L. Chernyak, A. R. Zhitnitsky and V. G. Serbo, JETP Lett. 26, 594 (1977) [Pisma Zh. 

Eksp. Teor. Fiz. 26, 760 (1977)]. 
[8] A. P. Bakulev, S. V. Mikhailov and N. G. Stefanis, Phys. Rev. D 67, 074012 (2003). 
[9] A. P. Bakulev, S. V. Mikhailov and N. G. Stefanis, Phys. Lett. B 578, 91 (2004). 
[10] V. M. Braun and I. E. Filyanov, Z. Phys. C 44, 157 (1989) [Sov. J. Nucl. Phys. 50, 511. 1989 

YAFIA, 50, 818 (1989 YAFIA,50,818-830.1989)]. 
[11] A. P. Bakulev and A. V. Radyushkin, Phys. Lett. B 271, 223 (1991). 
[12] A. P. Bakulev and S. V. Mikhailov, Z. Phys. C 68, 451 (1995). 

[13] A. P. Bakulev, S. V. Mikhailov and N. G. Stefanis, Phys. Lett. B 508, 279 (2001) [Erratum- 

ibid. B 590, 309 (2004)]. 
[14] A. P. Bakulev and N. G. Stefanis, Nucl. Phys. B 721, 50 (2005). 

[15] L. Del Debbio, M. Di Pierro and A. Dougall, Nucl. Phys. Proc. Suppl. 119, 416 (2003). 
[16] S. Dalley and B. van de Sande, Phys. Rev. D 67, 114507 (2003). 

[17] V. Y. Petrov, M .V. Polyakov, R. Ruskov, C. Weiss and K. Goeke, Phys. Rev. D 59, 114018 
(1999). 

[18] A. E. Dorokhov, JETP Lett. 77, 63 (2003) [Pisma Zh. Eksp. Teor. Fiz. 77, 68 (2003)]. 
[19] A. E. Dorokhov, Nuovo Cim. A 109, 391 (1996). 

[20] A. E. Dorokhov, M. K. Volkov and V. L. Yudichev, Phys. Atom. Nucl. 66, 941 (2003) [Yad. 
Fiz. 66, 973 (2003)]. 

[21] I. V. Anikin, A. E. Dorokhov and L. Tomio, gamma Phys. Lett. B 475, 361 (2000). 
[22] A. E. Dorokhov, M. K. Volkov and V. L. Yudichev, JINR-E4-2001-162. 

[23] S. V. Esaibegian and S. N. Tamarian, Sov. J. Nucl. Phys. 51, 310 (1990) [Yad. Fiz. 51, 485 
(1990)]. 

[24] M. Praszalowicz and A. Rostworowski, Phys. Rev. D 64, 074003 (2001). 

[25] M. Praszalowicz and A. Rostworowski, Phys. Rev. D 66, 054002 (2002). 

[26] E. Ruiz Arriola and W. Broniowski, Phys. Rev. D 66, 094016 (2002). 

[27] J. W. Chen, H. M. Tsai and K. C. Weng, Phys. Rev. D 73, 054010 (2006). 

[28] A. P. Bakulev, S. V. Mikhailov and N. G. Stefanis, Phys. Rev. D 73, 056002 (2006). 

[29] V. Y. Petrov and P. V. Pobylitsa, |hep-ph/9712203l 

[30] P. Ball and M. Boglione, Phys. Rev. D 68, 094006 (2003). 

[31] V. M. Braun and A. Lenz, Phys. Rev. D 70, 074020 (2004). 

[32] A. Khodjamirian, T. Mannel and M. Melcher, Phys. Rev. D 70, 094002 (2004). 

[33] P. Ball and A. N. Talbot, JHEP 0506, 063 (2005). 

[34] P. Ball and R. Zwicky, Phys. Lett. B 633, 289 (2006). 

[35] V. L. Chernyak, A. R. Zhitnitsky and I. R. Zhitnitsky, Nucl. Phys. B 204, 477 (1982) [Erratum- 

ibid. B 214, 547 (1983)]. 
[36] V. L. Chernyak, A. R. Zhitnitsky and I. R. Zhitnitsky, Sov. J. Nucl. Phys. 38, 775 (1983) 

[Yad. Fiz. 38, 1277 (1983)]. 
[37] E. V. Shuryak, Nucl. Phys. B 203, 93 (1982). 
[38] D. Diakonov and V. Y. Petrov, Nucl. Phys. B 245, 259 (1984). 
[39] D. Diakonov and V. Y. Petrov, Nucl. Phys. B 272, 457 (1986). 

[40] M. C. Chu, J. M. Grandy, S. Huang and J. W. Negele, Phys. Rev. D 49, 6039 (1994). 
[41] J. W. Negele, Nucl. Phys. Proc. Suppl. 73,92 (1999). 
[42] T. DeGrand, Phys. Rev. D 64, 094508 (2001). 

[43] P. Faccioli and T. A. DeGrand, Phys. Rev. Lett. 91, 182001 (2003). 

[44] P. O. Bowman, U. M. Heller, D. B. Leinweber, A. G. Williams and J. b. Zhang, Nucl. Phys. 



16 



Proc. Suppl. 128, 23 (2004). 

M. Musakhanov, Eur. Phys. J. C 9, 235 (1999). 

M. Musakhanov, |hep-ph/010"4T63| 

M. Musakhanov, Nucl. Phys. A 699, 340 (2002). 

P. V. Pobylitsa, Phys. Lett. B 226, 387 (1989). 

J. Gasser and H. Leutwyler, Nucl. Phys. B 250, 465 (1985). 

H. -Ch. Kim, M. M. Musakhanov and M. Siddikov, Phys. Lett. B 633, 701 (2006). 
H. -Ch. Kim, M. M. Musakhanov and M. Siddikov, in preparation. 
M. M. Musakhanov and H. -Ch. Kim, Phys. Lett. B 572, 181 (2003). 
S. i. Nam and H. -Ch. Kim, in preparation. 

N. G. Stefanis, W. Schroers and H. -Ch. Kim, Phys. Lett. B 449, 299 (1999). 
N. G. Stefanis, W. Schroers and H. -Ch. Kim, Eur. Phys. J. C 18, 137 (2000). 
A. Szczepaniak, C. R. Ji and S. R. Cotanch, Phys. Rev. D 49, 3466 (1994). 

M. Gockeler et al, |hep-lat/0510089j 

P. Ball, V. M. Braun and A. Lenz, |hep^pn /0603063 

V. M. Braun and I. E. Filyanov, Z. Phys. C 48, 239 (1990) [Sov. J. Nucl. Phys. 52, 126 (1990 
YAFIA,52, 199-213. 1990)]. 

P. Ball and R. Zwicky, JHEP 0602, 034 (2006). 

S. A. Gottlieb and A. S. Kronfeld, Phys. Rev. D 33 (1986) 227. 

G. Martinelli and C. T. Sachrajda, Phys. Lett. B 190, 151 (1987). 

D. Daniel, R. Gupta and D. G. Richards, Phys. Rev. D 43, 3715 (1991). 

A. R. Zhitnitsky, Phys. Lett. B 329, 493 (1994). 

S. i. Nam and H. -Ch. Kim, |arXiv:hep-p h/060504l"| 

V. L. Chernyak and A. R. Zhitnitsky, Nucl. Phys. B 201, 492 (1982) [Erratum-ibid. B 214, 
547 (1983)]. 



17 



